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PREFACE 


With  diminishing  world  fuel  supplies,  and  a  global  increase  in  fuel  price  over  the  last  ten  years,  the  reduction  of  aircraft 
drag  has  become  a  technology  of  major  importance  to  aircraft  manufacturers.  Likewise,  advances  in  test  and  evaluation 
techniques  have  facilitated  the  accurate  evaluation  of  drag  and  led  to  concurrent  developments  in  drag  prediction  methods. 
A  noteworthy  development  is  the  use  of  a  number  of  novel  Sow  control  methods  which,  through  either  passive  or  active 
interaction  with  the  flow  physics,  can  lead  to  substantial  drag  reductions. 

This  special  course  covers  some  of  the  more  recent  progress  in  drag  reduction,  measurement  and  prediction.  The  topics 
presented  discuss  the  different  sources  and  contributions  to  aircraft  drag  with  particular  emphasis  on  those  areas  in  which 
significant  new  developments  have  taken  place. 

The  course  begins  with  a  general  review  of  drag  reduction  technology.  Then  the  possibility  of  reduction  of  skin  friction 
through  control  of  laminar  Bow  is  discussed,  with  design  aspects  of  laminar  flow  control  hardware  included.  The  other 
possibility  of  skin  friction  reduction  through  modification  of  the  structure  of  the  turbulence  in  the  boundary  layer  is  also 
discussed. 

Methods  for  predicting  and  reducing  the  drag  of  external  stores,  of  nacelles,  of  fuselage  protuberances,  and  of  fuselage 
afterbodies  are  then  presented. 

Transonic  drag  rise,  the  prediction  of  viscous  and  wave  drag  by  a  method  matching  inviscid  flow  calculations  and 
boundary  layer  integral  calculations,  and  the  reduction  of  transonic  drag  through  boundary  layer  control  are  also  discussed. 
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COHPUT ATI ORAL  DRAG  ANALYSES  AND  MINIMIZATION; 
MISSION  IMPOSSIBLE? 
by 

J.W.  Slooff 

National  Aerospace  Laboratory  (NLR) 
Amsterdam,  The  Rather lands 


A  discussion  Is  given  of  the  possibilities  and  problems  associated  with  drag  prediction,  analysis  and 
minimization  by  naans  of  computational  fluid  dynamics  (CFD)  techniques  based  on  Invlscld  and  viscous- 
lnvlscld  Interaction  flow  models. 

Attention  Is  paid  to 

-  classical  and  computational  means  of  distinguishing  between  various  drag  components 

-  how  to  obtain  values  of  drag  components  from  currently  available  CFD  codes  and  the  problem  of  accuracy 

-  examples  of  computational  drag  diagnostics 

-  computational  drag  minimization 

The  discussion  Is  limited,  in  essence,  to  sub/transonlc  transport  (wing)  aerodynamics,  but  should,  at 
least  partly,  also  be  applicable  to  sub/transonlc  aerodynamics  of  fighter  aircraft. 

1  INTRODUCTION 

^  The  accurate  prediction  of  aircraft  aerodynamic  drag  Is  a  generally  recognized  and  respected  problem. 

So  la  the  accurate  measurement  of  drag  In  the  wind  tunnel  that,  eventually,  forms  the  basis  for  full  scale 
drag  prediction. 

In  the  past  15  years  Computational  Fluid  Dynamics  (CFD)  has  emerged  as  an  additional  and  complementary 
tool  for  aerodynamic  design  and  analysis.  The  purpose  of  this  lecture  Is  to  review  and  comsmnt  on  Its  role 
as  a  drag  prediction  and  analysis  tool. 

The  aerodynamic  design  process  of  aircraft  Is  characterized  by  a  sequence  of  design  and  anmlyaia 
cycles.  In  each  cycle  a  (further)  reduction  of  drag  will,  generally,  be  one,  but  not  the  only  objective. 
Identification  of  the  source  of  an  unacceptably  or  undesirably  high  drag  level  or  drag  variation  with  lift 
or  Mach  number  Is  a  prerequisite  for  a  successful  drag  reduction  program. 

Identification  of  drag  sources  may  follow  different  approaches.  The  classical  or  phenomenological  one 
la  based  on  the  availability  of  overall  force  (wind  tunnel)  data  only.  In  combination  with  simple,  semi- 
empirical  theory.  CFD,  as  we  shall  see  later,  offers  possibilities  for  a  more  physically/analytically 
oriented  approach  in  which  the  various  contributions  to  drag  are  distinguished  by  the  underlying  physical 
mechanisms  rather  than  by  the  observed  aerodynamic  force  variation  phenomena. 

The  main  part  of  this  lecture  Is  directed  towards  the  question  of  hov^th*7 various  components  of  drag, 
that  i  snjtf  4tftrt|5nS3,  in  imnli1 1 1  ail  liy,  mill  i  sn  ns  nrtr  tUned  ffffll  f  lie  imiunlly  available  CFD  codes.  As 
wlUrtMT^llus traced  through  examples,  accuracy.  In  this  context,  is  a  crucial  Issue. 

will  also  be  demonstrated  that,  in  spite  of  Its  current  shortcomings,  CFD  Is  a  powerful  tool  for 
drag  diagnostics.  The  final  part  of  the  lecture  contains  a  discussion  on  computational  drag  minimization. 

r 

2  DRAG  ANALYSIS,  INTRODUCTORY  REMARKS 

There  are  various  possible  ways  of  distinguishing  between  different  components  of  aerodynamic  drag, 

e.g. 

-  according  to  the  type  of  force  acting  on  the  body,  l.e.  preasure  or  friction  drag  (Fig.  1) 

-  according  to  the  physical  mechanisms  that  are  responsible  for  drag,  l.e.  friction,  causing  "viscous" 
drag,  vortex  shedding,  responsible  for  Induced  drag,  and  shock  wave  formation,  causing  wave  drag 
(Fig.  2). 

Viscous  drag  manifests  Itself  partly  as  friction  drag  and  partly  as  pressure  drag.  The  other  two  as 
pressure  drag  only.  Neither  of  the  two  classifications  given  above  Is  fully  manageable  in  a  situation  in 
which  only  total  force  data  are  available,  such  as  in  routine  (balance)  measurements  In  a  wind  tunnel.  For 
this  reason  the  classical  approach  to  drag  analyses  Is  based  on  the  observation  of  characteristic  features 
In  the  variation  of  drag  with,  e.g.,  lift  and  Mach  number.  We  may  call  this  the  phenomenological  approach. 
In  the  phenomenological  approach  distinction  la  usually  made  between 

-  (subsonic)  minimum  drag;  C_ 

mln 

-  (subsonic)  drag-due-to-llft;  CL 

L 

-  compressibility  drag;  AC^ 

compr 

-  trim  drag;  ACL 

trim 

-  miscellaneous  drag  (due  to  excressenses ,  etc);  ACQ 

misc 

In  other  words,  the  total  drag  Is  written  as 

<L  ■  C_  +  CL  +  AC_  +  AC_  +  ACL  (1) 

Dtot  min  L  compr  trim  wise 

or  a  similar  expression. 

The  (subsonic)  minimum  drag  (Fig.  3)  and  drag-due-to-llft  (Fig.  4)  are  usually  defined  at  some  suitable 
eubcrltlcal  free  stream  Mach  number  such  as  M  -  0.5  or  0.6.  Compressibility  drag  Is  defined  on  the  same 
basis  (Fig.  5). 

As  defined  above,  contains  moot  of  the  "viscous  drag"  and  some  Induced  drag.  However,  a  dletlnc- 
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tion  between  the  two  can  usually  not  be  Bade  unless  detailed  knowledge  about  the  span  loading  of  the  wing 
Is  available.  In  prelinlnary  design  studies  (L.  Is  usually  estimated  by  means  of  the  classical  "form 

■In 

factor"  methods  that,  e.g.  can  be  found  In  Boomer's  Book  "Fluid  Dynamic  Drag",  USAF  DATCOM  and  the  ESDU 

data  sheets. 

The  (subsonic)  drag  due  to  lift  contains  most  If  not  all  of  the  induced  drag,  some  viscous  drag  and, 
possibly,  some  wave  drag  at  the  higher  lift  values.  It  is  usually  estimated  or  fitted  by  means  of  the 
classical  induced  drag  formula 


(2) 


the  "Oswald"  factor  "e"  accounting  for  the  deviation  of  the  spanwlse  load  distribution  from  the  elliptical 
one  and  for  the  variation  of  the  viscous  drag  with  lift. 

Compressibility  drag  (Fig.  5)  contains  most  If  not  all  of  the  wave  drag,  the  Increase  of  viscous  drag 
with  Mach  number  and  possibly  some  Induced  drag  due  to  changes  in  span  loading  with  Mach  number.  A  general 
theoretical  basis  (apart  from  CFD  means)  for  the  prediction  of  wave  drag  does  not  exist. 

Trim  drag  contains  the  induced  drag  Increment  of  wing  plus  tail  due  to  tall  load  plus  some  viscous 
drag  and  wave  drag.  The  latter  two  because  the  wing,  at  least  for  conventional  configurations,  has  to 
operate  at  a  higher  lift  coefficient  than  untrlsmwd  due  to  the  negative  tall  load. 

The  classical  or  phenomenological  approach,  as  mentioned,  leade  to  a  manageable  way  of  distinguishing 
between  drag  components  In  case  only  overall  force  date  ere  available.  A  disadvantage  la,  however,  that  it 
often  remains  difficult.  If  not  lmposslbla,  to  point-out  the  precise  source  of  an  undesirable  dreg 
"surprise".  The  reason  being  that  these  classically  distinguished  drag  components  do  not  reflect  one  single 
but.  In  fact,  several  physical  mechanisms.  In  this  respect  an  analytical  approach  to  dreg  analysis,  that 
distinguishes  between  the  three  major  physical  mechanisms,  l.e.  friction  (viscosity),  vortex  shedding  and 
shock  wave  formation,  has  certain  advantages.  As  we  shall  see,  CFD  methods  may  serva  such  purpose. 


3  CLASSIFICATION  OF  CFD  METHODS  WITH  RESPECT  TO  DRAG  PREDICTION  POTENTIAL 


Prior  to  going  Into  the  details  of  the  possibilities  and  limitations  of  computational  drag  determina¬ 
tion  It  is  useful  to  sumarlze  briefly  which  drag  phenomena  ara  modelled  by  the  methematlcal/phyelcel 
models  underlying  the  currently  available  CFD  codes.  Flge.  6  end  7  may  serve  thle  purpose. 

Figure  6  presents  the  hierarchy  of  mathematical /physical  flow  modela  of  Interest  for  applied  CFD 
purposes.  Reynolds  averaged  Navler-Stokes  methods  can  provide  both  friction  drag  and  preaaura  drag  and,  at 
least  In  principle,  contain  all  Information  to  distinguish  between  viscous  drag.  Induced  dreg  sod  wave 
drag.  However,  because  of  the  very  large  computational  effort  required,  such  methods  are  not  (yat)  general¬ 
ly  available  In  the  applied  computational  environment,  at  least  for  problems  fnvolvlng  three  dimensions. 

For  this  reason  we  will  not  pay  any  further  attention  to  Nevler-Stokea  methods  in  this  lecture  note. 

On  the  next  level  (zonal  modelling)  we  must  distinguish  between  Invlscld  flow  end  viscous  flow 
(boundary  layer)  methods.  The  mutual  interaction  between  the  invlscld  and  viscous  parts  of  the  flow  field 
Is  modelled  through  an  Interaction  law  describing  tha  coupling  batween  the  two  zonae  at  a  suitably  chosen 
boundary. 

The  boundary  layer  methods  model  viscous  drag  only,  subject  to  the  boundary  conditions  provided  by  the 
external  invlscld  flow  through  the  Interaction  law. 

Three  further  levels  of  modelling  with  decreasing  physical  slgnlflcanca  can  ba  distinguished  for  lnvia- 
c Id  flow.  Methods  based  on  the  Euler  equations  model  both  Induced  dreg  end  wave  drag.  So  do  methods  based 
on  the  full  potential  equation  for  compressible  flow,  but  as  far  as  wave  drag  la  concamad.  In  an 
approximate  seaae  only.  The  approximation  Involved  Implies  a  limitation  to  weak  or  at  most  moderately 
strong  shocks  with  s  maximum  normal  Mach  number  of  about  1.3.  This,  however,  is  acceptable  for  many 
practical  purposes.  Figure  8,  taken  from  ref.  1,  provides  an  Impression  of  the  accuracy  of  the  potential 
flow  model  In  relation  to  tha  Eular  model  for  e  shock  Mach  number  of  about  1.3. 

Linear  potential  flow  methods,  such  as  psnal  methods ,  modal  Induced  drag  only  (plus  a  part  of  the  weve 
drag  In  case  of  linearized  supersonic  flow). 


4  DETERMINATION  OF  DRAG  FROM  CFD  CODES 


Generally  speaking  thara  are  two  different  waye  of  obtaining  drag  values  from  CFD  code  application. 

One  Is  the  direct  integration  of  presaura  and  friction  forces.  Me  will  call  this  the  near-field  approach. 
The  other  is  based  on  application  of  the  momentum  theorem  to  e  control  volume  enclosing  the  configuration. 
We  will  cell  this  the  far-fleld  approach. 

4. I  "Near-fleld"  assessment;  pressure  and  friction  drag 

The  accurate  determination  of  pressure  drag  from  wind  tunnel  pressure  measurements  is  a  well  raepected 
problem.  The  following  examples  illustrate  that  the  problsm  Is  of  comparable  magnitude  In  CFD;  at  least  In 
threa  dimensions,  and  If  ona  does  not  have  access  Co  a  sufficiently  fast  computer  with  sufficiently  large 

Figure  9  presents  an  example  of  the  pressure  drag  dependence  on  mash  d  nslty  for  a  2D  lifting  airfoil 
with  a  strong  shock  wave  as  calculated  with  a  "state  of  the  art"  Euler  code  employing  an  O-type  mesh  (N 
being  the  number  of  meeh  points).  The  results  indicate  that,  for  the  particular  case. considered,  a  196  x  36 
?  mesh  suffices  to  obtain  a  pressure  drag  coefficient  value  accurate  to  1  "count"  (10  j.  (The  letter  value 

seems  to  be  widely  accepted  as  a  "target"  accuracy  for  drag  measurement  in  wind-tunnel  testing).  Similar 
I  results  for  s  shock-free  supercritical  airfoil  as  obtained  with  two  different  Euler  codes  (of  which  one 

:  with  several  differant  distribution  of  mesh  points)  and  a  full  potential  code  are  presented  In  figure  10. 

|  All  these  codes  use  0-type  meahee  but  with  different  distributions  of  mash  points. 

The  message  of  figure  10  Is  that  the  sensitivity  of  the  results  for  mash  density  depends  strongly  on 
I  the  distribution  of  mesh  points.  For  Euler  cods  1  this  sensitivity  Is  clearly  unacceptably  high  for  the 


i 

i 
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particular  type  of  meah  uaed.  For  Euler  code  2  It  la  seen  to  vary  strongly  with  the  distribution  of  mash 
points.  For  the  full  potential  coda  a  192x96  wash  is  required  for  a  l  count  accuracy. 

It  can  be  concluded  fro*  figures  9  and  10  that,  with  a  suitable  distribution  of  mesh  points,  a  200x50 
to  200x100  mesh  will  generally  be  sufficient  to  obtain  pressure  drag  values  for  2D  airfoils  that  are, 
numerically,  accurate  to  1  count.  For  the  currently  available  general  purpose  type  of  main  frame  computer 
this  Is  certainly  feasible. 

An  aspect  that  may  also  affect  the  accuracy  of  CFD  pressure  drag  prediction  is  mash  extent.  The  point 
is  Illustrated  by  figure  11  *).  Shown  is  the  depence  of  the  pressure  drag,  as  computed  with  a  finite  volusw 
type  Euler  code,  on  the  distance  from  the  airfoil  at  which  the  far  field  boundary  condition  (In  this  case 
free  stream  plus  asymptotic  far-field  vortex  term)  is  prescribed.  The  particular  case  requires  a  mesh  ex¬ 
tent  of  about  100  chords  for  a  1  count  accuracy  I  (This  value  would  be  reduced  If  a  more  accurate  asymptotic 
far  field  boundary  condition  Is  used). 

A  3D  example  of  pressure  drag  dependence  on  mesh  density  is  presented  in  figure  12.  The  case  considered 
Is  an  application  of  a  modified  FL022  code  (kef.  3)  to  a  transport  type  wing  with  transonic  flow.  About  400 
relaxation  Iterations  were  performed  on  each  fine  mesh  plus  successive  mesh  refinement.  The  figure  Illu¬ 
strates  that,  for  the  particular  C-type  mesh  utilized  In  these  calculations,  something  of  the  order  of  540  k 
mesh  points  are  required  for  a  (numerical)  accuracy  of  1  count  if  the  calculations  are  performed  for  a  given 
lift.  The  situation  Is  distressingly  worse  if  the  calculations  are  performed  for  a  given  angle  of  attack. 

Examples  such  as  that  of  figure  12  serve  to  Illustrate  the  point  that  the  accurate  determination  of 
lnvlscld  pressure  drag  for  3D  wings  is  currently  (almost)  impossible  In  an  engineering  environment.  The 
general  experience  Is  that  this  applies  to  finite  difference/flnlte  volume  type  field  methods  as  well  as 
"panel"  methods.  However,  as  we  shall  see  later,  this  does  not  necessarily  mean,  that  CFD,  even  now,  cannot 
play  a  useful  role  in  drag  prediction/analyses  studies. 

The  problem  of  the  accurate  prediction  of  the  friction  drag  by  means  of  boundary  layer  codes  is  of  a 
somewhat  different  nature.  Unlike  the  (current)  situation  with  lnvlscld  flow  methods,  numerical  accuracy 
(truncation  errors)  is  not  the  major  source  of  concern  In  boundary  layer  calculations,  the  computational 
requirements  being  relatively  modest.  In  boundary  layer  computations  the  "absolute"  accuracy  associated 
with  the  particular  choice  of  mathematical/physical  model  Is  of  greater  Importance.  The  most  Important 
factors  involved  are: 

-  type  of  method  ("Integral"  of  "field"  method) 

-  transition  and  turbulence  models 

-  type  of  Interaction  law  used  ("weak"  or  "strong") 

-  type  of  flow  (attached  or  separated) 

Examples  of  skin  friction  prediction,  for  low  speed  flows,  are  presented  in  figures  13,  14.  In  one 
case  (Eppler  airfoil)  the  agreement  between  measured  and  calculated  skin-friction  distribution  is  seen  to 
be  quite  good.  In  the  other  case  (upper  surface  of  the  main  airfoil  of  a  hlgh-lift  configuration)  there  is 
substantial  disagreement  between  theory  (4  different  methods)  and  experiment,  the  reason  of  which  Is  not 
clear.  These  examples  Illustrate  that  the  accurate  prediction  of  skin  friction  may  not  (yet)  be  taken  for 
granted . 

It  Is  the  author's  general  Impression  that  the  "absolute"  errors  involved  in  friction  drag  prediction 
nay  be  as  large  as  5  to  15  counts  for  attached  flows  and  possibly  larger  for  separated  flows.  This,  how¬ 
ever,  does  not  exclude  the  possibility  to  obtain  meaningful  predictions  of  friction  drag  for  similar  types 
of  flow  on  similar  configuration  on  an  increment /decrement  basis. 


"Far-field  assessment  (momentum  theorem) 


Expressions  for  the  drag  acting  on  a  body  can  also  be  derived  by  applying  the  momentum  theorem  to  a 
control  vo luma  enclosing  the  body.  In  the  following  we  will  do  so  first  for  the  case  of  purely  lnvlscld 
flow  and  subsequently  for  lnvlscld  flows  with  "viscous"  lntersctlons  and  boundary  layer /wake  flows  ("zonal 
modelling") . 


4.2.1  lnvlscld  flow,  general  expressions 


Consider  a  control  voluma  of  the  type  as  Indicated  (for  a  half  wing)  in  figure  15.  In  the  far-field 
(Sw)  the  volume  is  bounded  by  a  plane  normal  to  tha  x-axla  far  upstream  of  the  wing,  a  similar  plane  (far) 
downstream  and  planes  parellsl  to  the  x-z  and  x-y  planes,  respectively,  each  at  a  large  distance  from  the 
wing.  In  the  near-field  the  volume  is  bounded  by  the  body  (wing)  surface  (S#)  and  (both  sides  of)  the 
shock-wave  (8  .  .  )  and  vortax  sheet  (S  .  )  surfaces.  The  reasons  for  excluding  the  vortex  sheet  from  the 

control  voluminfs  in  the  fact  Chat  the  ffcter  must  be  bounded  by  a  single  closed  surface.  Hence  a  cut  must 
be  Introduced  between  S  and  S.  . .  The  vortex  sheet  is  s  natural,  and,  as  we  shell  see  later,  sometimes 
nscesaary  choice  for  such  s  cud.  The  reason  for  excluding  the  shock  surface  will  become  apparent  shortly. 

Application  of  the  coneervatlon  law  for  momentum  to  the  control  volume  of  figure  15  leads  to  the 
Integral  expression 


+  (pq.n)u]  dS 


Sm  8  .  .  4-  8  .  +  S,  . 

B  shock  wake  (•>) 


p  Is  the  (local)  static  pressure,  q  is  the  total  velocity  vector,  p  la  tha  density,  n  is  tha  unit  outward 


normal  to  the  control  voluma  (hence  pointing  inside  S_,  S  .  .  and  S  .  )  and  u  is  the  component  In  the 

free-strsam  (x-)  direction  of  5.  9  sn0CK  **** 


*)  Figs.  9,  10  and  11  have  bean  derived  from  data  contained  by  reference  2. 


Splitting-off  the  body  pressure  integral  from  eq.  (3)  leads  to  the  following  relation  for  the  pressure  drag 
D  -  Opn  d  S  ■ 

p  77  * 
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This  can  be  rearranged  to  read 
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denotes  a  jump  of  the  particular  quantity  across  the  particular  surface  of  discontinuity. 

Note  that,  in  the  case  of  zero  normal  mass  flux  on  the  body  and  vortex  sheet,  which  we  will  consider  first, 
both  I  and  the  second  part  of  I-  are  zero.  Note  also  that  eq.  (5)  is  equally  valid  in  "full  potential" 
flow  and  "Euler"  flow.  However,  the  various  terns  distinguished  in  eq.  (5)  have  different  aeanings  in  these 
two  flow  models. 

In  the  (tine-dependent)  Euler  codes  that  are  currently  available,  aonentum,  as  well  as  aass  and 
energy,  are  conserved  in  the  whole  of  the  flow  field,  including  shocks  and  vortex  sheets.  The  latter  two 
are  captured  as  narrow  regions  with  high  velocity  gradient  rather  than  discontinuities  as  a  result  of  the 
miseries 1  and  artificial  dissipation  associated  with  the  discretization  schemes  that  are  being  used.  Within 
these  layers  nonentun  is  conserved  to  the  (numerical)  accuracy  implied  by  these  schemes.  It  follows  that  in 
the  Euler  flow  model  the  Integrals  I  and  I.  of  eq.  (5),  assuming  vanishing  truncation  errors,  are  identi¬ 
cal  to  zero.  As  already  mentioned,  the  fourth  integral  (I.)  is  also  zero  in  case  of  zero  normal  velocity  on 
SR.  Hence,  the  remaining  term  Ij,  i.e.  in  the  far-field  Integral,  contains  both  the  wave  drag  and  the 
Induced  drag. 

In  (conservative)  potential  flow  models  mass  as  well  as  energy  (and  entropy)  are  conserved  across 
shock  waves,  (or  rather  the  shock  layers  that  result  from  the  numerical  schemes),  but  momentum  normal  to 
the  shock  is  not.  The  latter,  in  fact,  is  the  very  mechanism  for  the  production  of  wave  drag  in  potential 
flow  models,  as  described  by  Stager  and  Baldwin,  Ref.  6  (see  also  Van  der  Vooren  and  Slooff,  Ref.  7).  It 
also  constitutes  the  reason  for  excluding  the  shock  surface  from  the  momentum  control  volume.  In  spite  of 
Che  non-physical  nature  of  this  seemingly  sC range  mechanism,  potential  flow  methods  generate  useful  numbers 
of  "wave  drag"  as  long  as  the  component  of  the  local  Mach  number  normal  to  the  shock  does  not  exceed  a 
value  of  say,  1.3.  The  point  was  already  illustrated  by  figure  8. 

From  the  discussion  Just  given  it  follows,  that  in  the  potential  flow  model,  the  Integral  Ij  in 
eq.  (5)  is  different  from  zero  and  represents  (an  approximation  to)  the  wave  drag.  The  far-field  Integral 
1^  now  contains  the  Induced  drag  only.  The  integral  I.,  as  mentioned,  is  zero  in  case  of  zero  normal  mass 
flux  on  S_.  Remains  to  discuss  the  "wake"  integral  I-* for  the  potential  flow  case. 

In  the  currently  available  potential  flow  codeszit  is  general  practice  to  model  the  vortex  sheet  as  a 
(curved)  plane  of  discontinuity,  coinciding  with  a  plane  of  Che  computational  mesh,  that  carries  a  Jump  in 
potential.  The  jump  In  potential  is  taken  to  be  constant  along  the  mesh  lines  in  this  plane  which  coincide, 
approximately,  with  the  free  stream  direction.  In  addition  the  difference  formulae  utilised  on  both  sides 
of  this  "vortex  sheet"  are  chosen  such  that  there  is  continuity  of  normal  mass  flux  across  the  "vortex 
sheet".  An  implication  of  this  is,  that  the  "flux"  term  of  momentum  is  continuous  across  the  vortex  sheet 
discontinuity  but  the  pressure  is  continuous  in  a  small  perturbation  sense  only.  Hence  the  term  I2  in 
eq.  (5),  although  small,  is  net  necessarily  identical  to  zero  in  potential  flow  calculations  with  a  priori 
selected  fixed  position  of  the  vortex  sheet.  The  author  is  not  aware  of  any  numerical  assessment  of  this 
term,  which,  indeed,  is  generally  neglected. 


4.2.2  The  evaluation  of  wave  drag  and  Induced  drag  from  potential  flow  codes 


In  the  preceedlng  paragraph  we  have  seen  that,  in  potential  flow,  the  wave  drag  can  be  determined,  at 
least  in  principle,  by  integrating  the  jump  in  momentum  across  the  shoe  surface (s) 7  (term  I.  of  eq.  (5)). 
In  practice  it  appears  that  the  accurate  determination  of  the  wave  drag  along  these  lines  poses  several 
difficulties. 

First  of  all,  as  noted  before,  shocks  are  smeared  by  numerical  capturing  schemes.  This  problem  can  be 
overcome  by  enlarging  the  contour  of  integration  to  ensure  that  the  shock  is  completely  enveloped  and 
calculating  the  loss  of  momentum  within  this  contour.  Figure  16  shows  permissible  contours  of  integration. 
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Th*  major  problem  associate d  with  calculating  the  wave  drag  along  tbesa  lines  la  tha  accurate  numeri¬ 
cal  evaluation  of  tha  Integral.  It  has  been  found  (Yu,  et.al..  Ref.  8)  that  very  fine  meshes  are  required 

to  perform  tha  integration  with  acceptable  accuracy  and  that  the  relaxation  process  for  solving  the  flow 
equations  must  be  continued  to  iteration  levels  that  are  far  beyond  those  of  current  engineering  practice. 
The  reason  for  the  latter  Is  In  the  fact  that  any  mass  flux  residual  (source  of  sink  like)  acts  as  a  spuri¬ 
ous  drag  of  thrust  source.  This  must  either  be  avoided  by  a  large  number  of  Iterations  or  be  corrected  for, 

see  reference  8. 

Making  use  of  the  mass  conservation  lav  applied  to  the  contour  S^ock  (Fig.  16)*),  l.e. 

-  $  p  5.5  d  s  -  »>hock  (6) 

Sshock 


where  R  .  .  represents  the  sum  of  the 

eq.  (5)  caS  oe  written  as 


residual  sources  in  the  flow  field  within  S  .  ,  , 

shock 
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the  Integral  Ij  of 
(7) 


With  the  second  term  of  (7)  representing  the  drag  or  thrust  associated  with  any  residual  sources  within 
^ahock*  first  term  represents  the  wave  drag,  l.e,: 
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[P  n  +  (u-u^)  p  q.n]  d  S  -  I 
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shock 


The  problem  of  accuracy  is  Illustrated  by  figure  17  which  was  taken  from  reference  8.  Por  a  transport 
type  of  wing-body  configuration  with  a  mesh  which  Is  typical  for  current  engineering  practice,  the  calcula¬ 
ted  wave  drag  varies  between  34  and  59  drag  counts,  depending  on  the  integration  contour  selected.  Evident¬ 
ly  much  finer  grids  (possibly  of  the  order  of  (.5  to  1)  x  10  grldpolnts)  are  required  to  calculate  the 
wave  drag  in  this  way  with  acceptable  accuracy. 

An  Interesting  alternative  way  of  evaluating  the  wave  drag  from  potential  flow  codes  has  been  proposed 
by  Garabedlan  et.al..  Ref.  9.  They  have  shown  that  the  Integration  of  the  artificial  viscosity  term  over 
the  supersonic  part  of  the  flow  field  provides  a  measure  of  the  wave  drag. 


Aa  discussed  In  section  4.2.1  the  Induced  drag  in  potential  flow  Is  represented  by  the  term  I3  of 
eq.  (5).  In  analogy  with  the  shock  Integral  Ij,  the  term  can  be  rewritten  as 

Ij  ■  -  ||  [p  +  (u-ua)  p  q.n]  d  S  +  u^  [  p  q.n  d  S  -  R^J  (9) 

$  s_+s  , 

(•)  B  wake 


with  defined  by 


(p  q.n)  d  S  +  p  q.n  d  S  •  R^ 
SB+Swake 

>nd  term  of  eq.  (9)  again  rep re 
' ,  the  induced  drag  is  given  by 

[p  nx  +  (u-u^)  p  q.n]  d  S  -  I* 


(10) 


With  the  second  term  of  eq.  (9)  again  representing  the  drag  acting  on  intentional  and/or  spurious  sources, 
(within  S^),  the  induced  drag  is  given  by 

(11) 


In  finite  volume  type  of  potential  flow  codes  this  integral  can  be  evaluated  on  suitably  chosen  grid 
boundary  planes.  The  experience  appears  to  be  (Ref.  8)  that  this  can  be  done  with  acceptable  numerical 
accuracy. 

An  alternative  la  to  let  the  boundary  surface  S. .  of  eq.  (11)  go  to  infinity  in  which  case  the 
integral  reduces  to  the  "Tref f ts-plane”  (x-+*)  Integral: 

Oj  ■  -  ||  [(p-p.)  +  (u-uj  pu]  d  S  (12) 

x-f» 


Assuming  small  perturbations  In  the  "Treff tz-plane", 
result  (ses  figure  18) 


Di-  |  1*1 


S  .  Ox 

wake 


wake 


d  t 


eq.  (12)  may  be  further  simp lifted  to  the  classical 


03) 


Por  details  of  the  derivation  and  approximations  Involved  see  e.g.  Lock,  Ref.  10.  Note  that  eq.  (13)  can  be 
ussd  to  calculate  the  drag  from  finite  dlfference/flnlte  volume  type  of  calculations  as  veil  as  from 
linear,  panel  type  methods. 

While  the  email  perturbation  assumptions  on  vhlch  eq.  (13)  is  based  are  violated  for  high  life  and  low 
aspect  ratio,  and  poealbly  at  the  (tip)  edge  of  the  vortex  sheet  in  any  case,  the  general  experience 
(Including  the  author *s)  la,  that  it  leads  to  Induced  drag  values  of  acceptable  accuracy  in  many  if  not 
moat  practical  cases. 


M  Mots  the  positive  direction  of  the  normals  on  sBkoch  and  S^t 
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4.2.3  Thg  evaluation  of  drag  from  Euler  codes 

Aa  discussed  in  section  4.2.1  the  Integral  I.  of  eq.  (5),  or  eq.  (9),  contains  both  the  wave  drag  and 
the  Induced  drag  In  "Euler  flow",  (as  veil  as  the  drag  associated  vlth  any  spurious  or  real  mass  sources 
inside  the  flov  field).  Reportedly  (Ref.  8),  the  numerical  evaluation  of  the  far  field  Integral  in  Euler 
codes  can  be  performed  with  acceptable  accuracy  In  2D,  but  not  In  3D.  The  main  problesis  seem  to  be  In¬ 
appropriate  (far-field)  boundary  conditions,  at  least  for  this  purpose,  and  lack  of  resolution  due  to  the 
relatively  small  number  of  mesh  points  that  can  be  afforded  In  engineering  applications. 

Splitting  the  total  Inviscid  drag  from  the  far-fleld  integral  into  wave  drag  and  Induced  drag  parts  Is 
not  a  trivial  matter.  As  discussed  by  Lock  (Ref.  10)  It  can  be  done  only  under  certain  simplifying  assump¬ 
tions.  To  the  author's  knowledge  it  has  not  (yet)  been  attempted  in  actual  3D  Euler  applications,  tn  the 
latter  the  problem  Is  further  exhanced  by  the  fact  that  vortex  sheets  are  captured  as  shear  layers  rsther 
than  discontinuities.  Further  exploration  of  the  possibilities  along  these  lines  deserves  attention. 


A  separate  evaluation  of  the  wave  drag  la  possible  by  considering  the  entropy  production  of  shock 
waves,  or,  rather,  the  entropy  produced  within  a  contour  enveloping  the  numerically  smeared  shock  layers 
(Fig.  19). 

For  this  purpose  use  can  be  made  of  an  expression  originally  derived  by  Osvstitsch  (Ref.  II) 


D 

wave 


Tw 
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(S-S^)  p  q.n  d  s 


Sshock 


(14) 


(S  denoting  entropy  and  Tw  the  freestresm  temperature),  or  variants  thereof  (see  also  Ref.  10). 

The  experience  appears  to  be  (Ref.  8)  that  eq.  (14)  can  be  evaluated  with  acceptable  accuracy  In  2D 
but  not  (yet)  In  3D.  The  main  problems  being  spurious  entropy  production  (dissipation)  associated  with  the 
numerical  schemes.  In  particular  near  the  flov  boundaries,  and  lack  of  resolution  (mesh  points). 

Sumarlzlng  the  discussion  on  Euler  codes  we  may  conclude  that  the  problem  of  determining  drag  vlth 
acceptable  accuracy  Is  even  more  severe  than  for  potential  flov  codes.  In  particular  If  distinction  between 
wave  drag  and  Induced  drag  Is  the  objective. 


4.2.4  Zonal  modelling 

Assessing  the  various  components  of  drag  becomes  more  complicated  when  the  effects  of  viscosity  In  the 
flov  field  are  taken  Into  account  by  matching  the  outer  lnvlscld  part  of  the  flov  field  with  the  (inner) 
viscous  or  boundary  layer  part  through  an  appropriate  Interaction  law  (zonal  modelling).  It  appears,  see 
Llghthill,  Ref.  12,  that  there  are  several  alternative  ways  for  matching.  However,  the  moat  convenient,  and 
Indeed  moat  cosaaonly  adopted  possibility,  la  to  simulate  the  displacement  effect  of  the  boundary  layer  and 
wake  on  the  outer  lnvlscld  flow  through  Imposing  a  "transpiration  mass  flux"  boundary  condition  on  the  body 
surface  and  wake  centerline  or  vortex  sheet.  In  the  "weak"  interaction  formulation,  which  la  applicable  to 
attached  flows  only,  the  transpiration  mass  flux  at  the  body  surface  and  on  the  vortex  sheet  In  the  invis¬ 
cid  flow,  and  the  tangential  mass  flux  at  the  edge  of  the  boundary  layer  In  the  boundary  layer  flow,  are 
chosen  such  that  the  normal  and  tangential  mass  flux  components  in  both  the  lnvlscld  and  boundary  layer 
flow  match  at  the  location  S  ,  figure  20,  of  the  ed^  >f  the  boundary  layer  and  wake.  In  the  "strong" 
interaction  formulation,  whlcR*?a  (also)  applicable  to  separated  flows,  the  sequence  of  "inviacid"  and 
"viscous"  calculations  (interaction  law)  la  chosen  differently  for  reasons  of  computational  stability  and 
efficiency,  with,  however,  the  same  objective  of  matching  the  inviscid  and  boundary  layer  flow  at  the  edge 
of  the  boundary  layer  (see,  e.g.  Refs.  13,  14). 

With  zonal  modelling  in  the  sense  described  above,  we  have,  in  the  lnvlscld  part  of  the  flow  problem, 
apart  from  the  actual  lnvlscld  flow  outside  the  surface  S  ,  an  "equivalent"  Inviacid  flow  inside  of 
S  .  .  In  the  real  v4*»coua  flov,  the  region  within  S  .  fa?*of  course,  occupied  by  the  boundary  layer  and 
wall.  In  case  the  boundary  layer  and  wake  flow  are  modelled  by  the  "classical"  boundary  layar  equations 
there  are  no  pressure  gradients  across  the  boundary  layer  and  wake.  However,  triple  deck  theory  for  turbu¬ 
lent  boundary  layers  (for  a  survey  see,  e.g.  Ref.  14)  teaches  that  normal  pressure  gradients  are  not  negli¬ 
gible  near  the  trailing  edge,  near  separation  lines  and  where  shock  waves  Intersect  the  boundary  layer.  The 
formulations  given  below  are  chosen  such  that  they  also  cover  the  case  with  normal  pressure  gradients  in 
the  viscous  flov. 


With  the  purpose  of  relating  "edge"  and  body /wake  quantities  we  will  now  apply  the  momentum  theorem 
first  to  the  equivalent  Inviscid  flov  within  S  (Fig.  20). 

It  reads  ,dge 

tp  +  (p  q.n)  u)  d  S  •  0  (15) 

S.+S  .  +S  .  +8S  .  ,+dS,  . 

B  wake  edge  shock  («*) 


Note  chat,  for  the  same  reasons  as  discussed 
the  area  of  integration. 

Distinguishing  the  various  parts  of  the  area 


in  sections  4,2.1  and  4.2.2 
of  integration  eq.  (15)  can 


we  have  excluded  the  shock (a) 
be  written  aa 


from 
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As  before,  we  will  consider  the  weaning  of  the  various  terms  in  eq.  (28),  In  particular  for  the  situation 
that  the  downstream  boundary  (6)S^  of  the  control  volume  tends  to  infinity. 

It  follows  froa  the  proceeding  discussions  that  in  case  the  outer  lnvlscld  flow  is  aodelled  by  the 

potential  equation  the  wave  drag  Is  aodelled  by  the  terae  l!-6l!  _  +  dl.  .  end  the  Induced  drag  by  the 

term  TTT  Considering  the  tera  I_  It  nay  be  noted  that  in  two-dimdnsional  rlow*§*u„  for  x+*%  so  that  I,*0. 

In  three  dlaenslona  with  lift  tnis  is  not  exactly  the  case  because  of  the  flow  perturbations  Induced  By  the 
trailing  vortices.  However,  as  Indicated,  a.o.  by  Lock  (Kef.  10),  the  aagnltude  of  the  factor  (u-uj ,  Is 
generally  negligible,  except,  perhaps  for  high  lift  and  saall  aspect  ratio.  It  is  further  noted,  that  the 
factor  (p  u-pU),  when  integrated  over  68.  . ,  is  of  the  order  of  the  dlsplaceaent  voluae  of  the  wake. 
Hence  I_,  H  all  probability,  can  be  safely ^neglected  in  aost  practical  cases. 

The  tera  1^,  obviously,  represents  the  vlacoua  drag.  We  will  return  to  this  tera  in  section  4.2.5. 

In  case  the  outer  lnvlscld  flow  is  aodelled  by  aesns  of  the  Euler  equations  we  have,  as  mentioned 
before,  1,-81.  ,  +  6l,  .  -  0,  while  I'  contains  both  the  Induced  drag  and  the  wave  drag  of  the  equi¬ 
valent  invlscid*?Yov.  Wnlie?  +0  for  potential  flow,  at  least  In  two  dlaenslona,  it  is,  in  general,  *  0  In 

case  of  Euler  flow  because  u  *  um  in  the  wake  of  shock  waves.  Hence,  I_,  like  the  tera(s)  61  in  potential 
flow,  can  be  considered  as  a  (viscous)  correction  to  the  lnvlscld  wave  <lrag.  It  takes  account  of  the  fact 
that  the  wave  drag  associated  with  those  parts  of  the  shock  waves  of  the  equivalent  lnvlscld  flow  within 
S  .  should  be  subtracted  from  the  total  "lnvlscld"  wave  drag.  The  author  Is  not  aware  of  any  quantitative 
o!u8$merical  assessment  of  this  term.  The  term  I.  again  models  the  viscous  drag. 

It  should  be  noted  that  if  the  downstream  end  of  the  control  volume  Is  taken  at  a  finite  distance  from 
the  configuration  it  Is  no  longer  possible  to  Identify  the  various  drag  sources  (Induced /vortex,  wave, 
viscous  drag)  directly  and  uniquely  with  the  terms  of  eq.  (28). 


4.2.5  Viscous  drag 

As  mentioned  above,  the  viscous  or  boundary  layer  drag  Is  represented  by  the  term  1^  of  eq.  (28),  or 

Dvl.c  *  11  (u_U)  0  U  <*  S  <-  V  (29) 


Evaluation  of  eq.  (29)  requires,  in  the  general  case,  a  method  for  the  calculation  of  three-dimensional 
wakes.  The  problem  Is,  however,  that,  even  for  planar,  l.e.  non-rol ling-up  wakes  such  methods  are  not 
generally  available. 

Limiting  the  discussions  to  wings,  the  approximate  procedure  that  Is  usually  followed  to  circumvent 
this  problem  Is  to  compute  the  boundary  layer,  and.  If  possible,  the  wake,  up  to  the  trailing  edge  or  some 
finite  distance  downstream  of  the  trailing  edge  and  to  extrapolate  the  development  of  the  viscous  wake  to 
downstream  infinity. 

For  this  purpose,  eq.  (29)  Is  replaced  by 

\i.c  - {  °.<T')  *r  (30) 

where 
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la  the  (local)  momentum  thickness  at  downstream  Infinity.  (y',z')  are  coordinates  along  and  normal  to  the 
wake  in  the  plane  x— +“.  Note  that  the  expressions  (30),  (31)  are  applicable  only  In  the  case  of  planar, 
l.e.  wlng-llke  wakes. 

In  two-dimensional  flow,  the  classical  Squire  and  Young  method  (Ref.  15)  or  variants  thereof  are 
usually  applied  to  relate  to  the  momentum  thickness  9(x)  calculated  by  means  of  a  boundary  layer/wake 
method  at  some  finite  distance  from  the  trailing  edge.  Under  the  assumption  that  the  9hape  function  H 
varies  linearly  with  In  U  ,  (U  being  the  velocity  at  the  edge  of  the  boundary  layer),  the  Squire  and  Young 
method  leads  to  the  expreSsion® (In  "boundary  layer"  notation) 


(32) 


The  procedure  used  for  3-D  wings  Is  generally  similar  to  the  one  described  above.  In  this  case  8(x)  Is 
usually  taken  as  8-g  at  the  trailing  edge  and  some  extended  form  of  the  Squire  and  Young  method,  derived 
for  Infinite  yawed  wings,  Is  used  at  each  span  station  to  relate  to  (see  e.g.  Lock,  Ref.  10). 

It  is  appropriate  at  this  point  to  stress  the  point  that  the  Basic  assumptions  underlying  the  Squire 
and  Young  method  are  not  valid  In  case  of  (highly)  asymmetric  wakes  (aft-loaded  airfoils),  wakes  in 
(strong)  pressure  gradients  and  highly  three-dimensional  wakes.  In  such  conditions  methods  of  the  type 
sketched  above  should  be  applied  with  caution. 

Figures  22,  23  illustrate  results  of  an  application  of  the  procedure  to  a  3D  transport  type  wing  and  a 
comparison  with  experimental  (wake  rake)  data  (Ref.  16).  A  fully  three-di  ensional  boundary  layer  code  was 
used  to  calculate  6  upto  the  trailing  edge. 

Figure  23  illustrates  that  even  for  a  relatively  simple  case  as  that  of  figure  22  the  difference  between 
calculated  and  measured  viscous  drag  may  be  substantial.  Although  there  are  numerous  potential  sources  of 
error  In,  at  least,  the  computational  methods  and  procedures  Involved,  the  precise  origin  of  such  discre¬ 
pancies  la  by  no  means  clear. 

An  example  of  viscous  drag  prediction  for  a  fairly  complicated  2D  case  (airfoil  with  flap  at  low 
speed,  from  Ref.  17)  Is  reproduced  In  figure  24.  The  figure  serves  the  purpose  of  Illustrating  that  turbu¬ 
lence  modelling  Is  of  crucial  Importance  for  wakes  In  strong  longitudinal  pressure  gradients  such  as  occur 
In  the  main  wake  of  flapped  airfoils. 
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4.3  Cross-checking  near-field  and  far-flald  values 

Is  view  of  the  Inaccuracies  end  uncertainties  discussed  and  illustrated  above,  it  is  advisable  to  cow- 
put  a  the  different  conponents  of  drag  In  as  many  ways  as  possible  and  compare  the  reaulta  (vis.  figure  25). 
Differences  between  drag  component  values  coaputed  In  different  ways  may  provide  an  indication  of  the  error 
levels  that  are  possibly  Involved. 

On  the  inviscld  side,  an  obvious  check  Is  to  cowpare  the  Integrated  (near-fleld)  pressure  drag  against 
the  sun  of  the  Induced  drag  and  wave  drag  as  obtained  from  the  (far-fleld)  momentum  considerations.  Eq.  (5) 
represents  the  basic  equation  for  such  a  check. 

In  case  of  potential  flow  It  nay  be  wore  convenient  to  rewrite  eq.  (5)  In  terms  of  the  expression  (7), 
(8)  and  (9)-(13)  for  the  (corrected)  wave  drag  and  Induced  drag. 

This  yields 
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with  D  given  by  eq.  (8)  and  D  by  eq.  (11),  (12)  or  (13).  The  third  term  on  the  right  hand  aide  of 
eq.  OJfjy®is  non-zero  In  case  of  surface  (and  wake)  "transpiration”  only.  It  can  be  Interpreted  as  the 
pressure  drag  (of  the  equivalent,  Inviscld  flow)  Induced  by  the  viscous  displacement  effect  (D  ) .  The 
fourth  term,  D  ,  represents  the  effect  of  a  pressure  jump  across  the  Inviscld  wake  discontinuity 

pwake 

(vortex  sheet).  It  is  zero  in  case  of  a  fully  relaxed  vortex  sheet  and  vanishing  viscous  wake  thickness  or 
viscous  wakes  with  zero  normal  pressure  gradients  (see  also  section  4.2.1).  The  terms  R,  .  and  R  .  ,  re¬ 
present  the  effect  of  possible  residual  sources  inside  the  flow  field  (section  4.2.2).  TCey  are ,o?Ccourse , 
zero  in  case  of  a  fully  converged  solution. 

In  case  the  Inviscld  flow  Is  modelled  by  the  Euler  equations  It  follows,  from  the  discussion  of 
sections  4.2.1  through  4.2.3,  that  the  corresponding  expression  Is 
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with  I’  (eq.  (11))  containing  both  the  wave  drag  and  the  Induced  drag. 

In  viscous  flow  a  similar  check  can,  of  course  be  performed  by  computing  the  total  drag  by  Integrating 
the  surface  pressure  and  skin  friction,  eq.  (20),  as  well  as  by  eq.  (28). 

An  additional  checking  possibility  Is  obtained  by  eliminating  I’  from  eqs.  (28)  and  (33),  (34), 
leading  to  the  general  expression. 
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In  case  of  potential  outer  flow  we  have  I.+0,  at  least  In  two  dimensional  flow  (cf.  section  4.2.4). 
In  case  of  "Euler”  outer  flow  this  reduces  to 
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Rote  that.  In  the  absence  of  shock  waves,  with  zero  normal  pressure  gradients  across  boundary  layer  and 
wake,  and  with  zero  mass  flux  residual,  eqs.  (37),  (38)  reduce  further  to 
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(39) 


expressing  Chet  Che  total  viscous  dreg  sinus  the  friction  dreg  equals  (should  equal),  the  pressure  drag  due 
to  viscous  effects. 


4.4  Englnearln 


p roach  to  computational  drag  analysis 


In  the  preceeding  sections  we  have  seen  that  we  are  still  far  from  a  situation  In  which  ve  can  compute 
drag  vlth  acceptable  absolute  accuracy.  An  additional  problem  Is  that  the  majority  of  the  available 
lnvlscld  flow  codes  provide  values  for  Integrated  pressure  drag  only,  so  that  a  division  Into  wave  drag  and 
Induced  drag  Is  not  readily  available. 

In  spite  of  this  situation  It  has  proven  to  be  possible  to  devise  engineering  procedures  for  deter- 
mining  drag  from  CFD  codes  with  usable  accuracy;  at  least  in  applications  where  drag  increments,  due  to 
(small)  changes  in  flow  conditions  or  geometry,  rather  than  absolute  accuracy  Is  the  subject  of  Interest. 
Examples  of  such  procedures  may  be  found  in  a.o.  Refs.  3  and  18. 

In  Ref.  3  the  Induced  drag  of  a  wing-body  configuration  Is  calculated  from  Tref ftz-plane  considera¬ 
tions,  utilizing  the  spanvlse  distribution  of  circulation  obtained  from  a  potential  flow  code  and  assuming 
that  there  Is  no  vortlclty  shed  from  the  body. 

The  following  procedure  Is  used  to  obtain  an  estimate  for  the  wave  drag.  First,  the  total  lnvlscld 
drag  of  the  wing  as  determined  from  Integrated  pressures  Is  compared  with  the  Induced  drag  of  the  wing  for 
a  number  of  subcritical  flow  conditions.  Assuming  that  the  Induced  drag  Is  determined  with  sufficient 
accuracy  this  comparison  yields  a  correction  AC^  for  the  pressure  drag  which  Is  expressed  as  a  function  of 

P 

wing  lift  and  Mach  number  as  follows, 
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The  quadratic  C  -term  In  (40)  Is  suggested  by  the  fsct  that  the  Induced  drag  (and  hence,  in  subcritical  flow, 
also  the  pressure  drag)  Is  a  quadratic  function  of  the  lift;  one  would  expect  the  error  In  this  to  be  propor¬ 
tional  to  the  level  of  the  induced  drag.  The  terms  with  M*  and  (1  -  M*)  Yare  suggested  by  the  Raylelgh- 
Janzen  M^-expanslon  theory  and  Prandtl-Glauert  theory,  respectively.  The  cross-term  M*c£  Is  merely 
empirical.  wing 


The  constants  C.  to  C,  are  determined  by  equating  the  corrected  pressure  drag  and  the  far-fleld 
induced  drag  for  six  subcritical  flow  conditions  involving  three  different  C  -values  and  three  Mach 
numbers.  The  experience  is  that  In  this  way  the  corrected  pressure  drag  In  the  whole  subcritical  flow 
regime  can  be  determined  with  an  accuracy  of  a  few  drag  counts.  With  the  constants  to  C  known  the 
correction  formula  (40)  Is  also  used  to  correct  the  pressure  drag  In  supercritical  flow  conditions.  For 
supercritical  flow  conditions  the  wave  drag  Is  then  determined  by  subtracting  the  Induced  drag  of  the  wing 
from  the  corrected  pressure  drag,  l.e. 


CD  "  S  +  (H-’  CL  4  > 

wavc  P„lno  P  wln8 


Note  that  the  procedure  assumes  that  wave  drag  Is  generated  by  the  wing  only. 

An  alternative  for  the  "pressure  drag  error  extrapolation"  technique  described  above  may  be  to  use  a 
mesh  density  extrapolation  procedure  of  the  type  illustrated  by  figure  12. 

The  pressure  drag  error  extrapolation  procedure  of  Ref.  3  forms  part  of  a  fairly  complex  system  of 
computer  programs  for  the  prediction  of  the  aerodynamic  characteristics  of  wing-body  configuration  In 
transonic  flow  (Fig.  26).  In  this  system  wing  viscous  drag  is  computed  by  a  fully  3D  finite  difference  type 
of  boundary  layer  code  in  combination  with  an  extended  form  of  the  Squire  and  Young  method  (section  4.2.5). 
Figure  27  provides  comparisons  of  measured  and  computed  drag  for  a  narrow-body  transport  configuration,  ob¬ 
tained  through  the  system  of  figure  26,  illustrating  the  usefulness  of  the  drag  prediction  procedure.  Note 
that  the  variation  of  drag  with  lift  and  Mach  number  is  predicted  reasonably  well  as  long  as  the  flow  Is 
attached;  even  the  absolute  drag  level  is  fairly  close  if  the  computed  wing  drag  Is  added  to  the  (viscous) 
drag  of  the  body  alone  as  measured  In  the  wind  tunnel.  Further  examples  may  be  found  In  ref.  3. 

The  pressure  drag  extrapolation  technique  may  also  be  used  wing  section  wise  to  obtain  estimates  for 
the  distribution  of  wave  drag  along  the  span  of  the  wing.  Additional  Information  provided  by  the  system  of 
figure  26  contains  boundary  layer  and  wing  section  viscous  drag  quantities  (Fig.  28)  and,  of  course, 
detailed  pressure  distributions. 


5.  COMPUTATIONAL  DRAG  MINIMIZATION 

The  past  decade  has  seen  repeated  efforts,  some  (partly)  succesful,  others  less  so,  to  directly 
address  the  problem  of  drag  minimization  by  combining  computer  codes  for  aerodynamic  (drag)  analyses  (flow 
solvers)  with  numerical  optimization  algorithms.  The  objective  being  that  the  optimization  algorithm  con¬ 
trols  variations  of  a  number  of  Independent  variables,  such  as  parameters  defining  the  geometry,  with  the 
purpose  of  finding  the  particular  combination  of  parameters  that,  subject  to  given  constraints,  leads  to  a 
minimal  drag.  In  this  process,  the  function  of  the  flow  solver  Is  to  provide  drag  values  for  each  combina¬ 
tion  of  values  of  the  Independent  variables  that  Is  considered  to  be  feaslb  e  and  "Interesting"  by  the 
optimization  algorithm. 

Generally  speaking  computational  (or  numerical)  drag  minimization  requires 

-  a  choice  of  object  function  (In  this  case  the  drag*  but  other  choices  are  possible) 

-  a  flow  solver  (aerodynamic  analysis  code) 

-  a  choice  of  Independent  variables 

-  a  choice  of  constraint  functions  defining  that  part  of  the  solution  space  that  Is  considered  to  be 
feasible  from  the  engineering  or  another  (e.g.  numerical)  point  of  view 
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-  ah  optimisation  algorithm. 

With  respect  to  the  choice  of  Independent  variables  one  may  distinguish  two  different  approaches .  One, 
and  Indeed  the  most  common  choice  is  to  use  parameters  defining  the  geometry  as  the  Independent  variables. 
This  requires  a  direct  or  analyses  type  of  flow  solver  only.  The  approach  is  generally  referred  to  as 
direct  numerical  optimisation. 

While  application  of  the  direct  numerical  optimisation  approach  have  seen  the  use  of  various  aero¬ 
dynamic  codaa,  the  feasible  directions/gradient  optimisation  algorithm  CONMDt/COPES,  developed  by 
Vanderplaata  (Ref.  19)  seems  to  be  used  elmoet  exclusively.  In  particular  In  combination  with  transonic 
flow  codaa. 

The  approach,  pioneered  by  Hicks  at  al  (taf.  20)  owe  Its  existence  entirely  to  the  availability  of 
large  and  feat  computer  systems,  because  of  the  excessively  large  computational  requirements ,  at  least  In 
3D,  the  approach  is  soma times  referred  to  as  "design  by  bruts  force".  Nevertheless  it  holds  great  potential 
for  the  future.  A  reappraisal  of  the  technique  has  bean  given  by  Hicks  (Ref.  21).  Several  chastening 
experiences  are  also  reported  In  Ref.  28. 

A  generalised  flow  diagram  of  the  numerical  optimisation  technique  Is  presented  in  figure  29.  The 
process  la  Initiated  by  the  choice  of  an  aerodynamic  object  function  F  that  la  to  be  minimised  (In  this 
case  the  drag)  a  number  of  quantities  to  be  constrained  G  end  e  eet  of  design  variables.  The  constraints 
can  be  of  aerodynamic  of  geometric  nature;  e.g.  C.  and/oirt/c  greater  than  a  specified  value.  The  design 
variables  are  generally  taken  to  be  the  coefficients  A^  of  s  number  of  shape  functions 

n 
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describing  (modifications  to)  the  starting  wing  geometry. 

The  process  begins  by  perturbing.  In  sequence,  each  of  the  shape  function  coefficients  A..  The  result¬ 
ing  n  shapes  are  analyzed  by  means  of  the  aerodynamic  program  (determination  of  F  and  G  *e)  and  the  derl- 
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vet Ives  5I~»  5.  •  01  rather  the  difference  quotients  tt~  are  determined.  The  next  step  is  the 
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formation,  by  the  optimization  program,  of  the  gradient  ^F  and  the  determination  of  the  direction  of 
steepeet  descent  of  F,  in  the  n-dlmenslonal  space  formed  by  the  basis  vectors  A  ,  while  satisfying  the 
constraints.  The  optimisation  program  than  executes  s  number  (typically  3)  of  steps  In  this  direction,  with 
another  aerodynamic  analysis  performed  at  each  step,  until  either  s  constraint  la  met  or  F  attains  s 
minimum.  In  the  first  esse,  or  when  the  minimum  of  F  Is  lower  than  the  previous  minimum,  the  process  Is 
repeated;  new  gradients  are  determined,  etc.  When  the  latest  minimum  or  F  Is  equal  to  or  higher  then  the 
previous  one  the  process  Is  terminated. 

The  optimization  process  described  above  requires  typically  10  complete  cycles  or.  In  other  words, 
10(n+3)  analysis  calculations  (Ref.  19).  This  immediately  illustrates  the  weakes-  point  of  the  numerical 
optimization  approach.  In  order  to  keep  the  computational  effort  required  within  reasonable  bounds  one  has 
to  put  severe  limitations  on  the  number  n  of  design  variables.  In  particular  In  3D  flow.  The  problem  is 
enhanced  by  the  feet  that  for  acceptable  convergence  of  the  optimization  process  it  Is  necessary  to  avoid 
"numerical  noise"  In  the  partial  derlvates  of  the  object  function  (Refs.  23,  24).  This  requires  that  the 
relaxation  process  in  each  analysis  calculation  must  be  continued  until  the  residual  has  reached  a  level 
beyond  that  which  la  often  customary  in  "normal"  analysis  calculations.  It  also  appears  to  exclude  the  use 
of  analysis  codes  with  simple  boundary  layer  corrections  (Ref.  21).  The  reason  for  the  latter  is  that  the 

airfoil  aerodynamic  quantities  do  not  vary  consistently  enough  when  boundary  layer  and  potential  flow  are 

coupled  In  the  weak  interaction  aenaa. 

One  way  to  reduce  the  number  of  analysis  calculations  required  in  3D  applications  la  to  evolve  the 
design  variables  In  a  series  of  steps  (Ref.  26).  For  example  by  first  designing  the  upper  surface*  section 
by  section,  going  from  root  to  tip  end  then  the  lower  surface.  Clearly  It  la  also  important  to  select  a 

starting  geometry  having  aerodynamic  characteristics  which  are  already  does  to  the  target.  This  asks  for 

an  Information  system/data  base  approach.  With  previous  experience  stored  in  the  data  base,  the  latter  can 
be  searched  for  the  moat  aultabla  starting  solution.  As  described  in  Ref.  19  the  data  bass  approach  can 
also  be  used  to  speed-up  the  convergence  of  numerical  optimization  by  at  least  a  factor  two.  With  the 
results  of  all  preceding  geometry  perturbations  stored  it  is  possible  to  construct  higher  partial  deriva¬ 
tives  of  the  object  function  sod  utilize  higher  order  gradient  methods. 

With  the  severe  limitations  on  n,  the  choice  of  the  shape  functions  is  of  utmost  importance.  The 
choice  should  be  directed  towards  describing  a  sufficiently  wide  class  of  practical  solutions.  While  simple 
polynomial  expressions  were  used  In  early  applications  (Refs.  20,  22)  of  the  numerical  optimization  con¬ 
cept,  a  more  sophisticated  class  of  ahape  functions  describing  mors  local  geometry  modifications  mas  used 
in  later  applications*  (Refs.  23,  24,  25),  However,  as  discussed  in  Ref*  26,  there  Is  a  need  for  still 
batter  shape  functions  with  even  more  localized  curvature  variations.  In  fact  it  can  be  argued  that  while 
curvature  based  shape  functions  are  suitable  for  arses  with  subcritlcal  flow,  slope  based  shape  functions 
might  be  more  appropriate  in  areas  with  locally  supersonic  flow. 

An  Interesting  choice  for  the  shape  functions  la  discussed  by  Aldale  at. si.  (Ref.  29).  They  consider 
shape  functions  generated  by  means  of  feeding  certain  pressure  distribution  modifications  into  an  inverse 
program.  The  result  Is  a  sat  of  design  shapes  that  are  (almost)  orthogonal  in  an  aerodynamic  aenee,  that 
la,  affect  only  one  specific  pressure  distribution  characteristic  and  no  other  ones. 

While  the  choice  of  the  design  variables  la  of  grant  practical  significance,  the  precise  choice  of  the 
object  function.  In  conjunction  with  the  choice  of  the  aerodynamic  and  geometric  constraints,  la  of  both 
more  fundamental  and  practical  Interest.  In  two-dimensional  transonic  applications  (Rafs.  19,  20*  22,  23) 

It  has  been  customary  to  minimise  the  wavs  drag  subject  to  constraints  on,  e.g.,  airfoil  thickness  or 
volums,  lift  and/or  pitching  moment.  Although  It  is  clear  that  constraints  are  necessary  In  a  meaningful 
drag  minimization  problem  It  la  by  no  means  clear  how  exactly  the  problem  should  be  formulated  In  order  to 
guarantse  a  unique  solution.  The  problem  is  Illustrated  by  figure  30,  taken  from  Ref.  19.  Shown  Are  the 
results  of  two  drag  minimisation  runs  with  Identical  free  stream  conditions  and  Identical  constraints  on 
lift  and  airfoil  volume.  Only  the  starting  solutions  differ.  As  illustrated  by  the  figure  the  two  resulting 
airfoils  are  totally  different  In  shape.  Clearly  the  problem,  as  formulated,  has  more  than  one,  local 
minimum  and  neither  of  the  two  necessarily  represents  the  absolute  minimum.  An  Interesting  discussion  on 
criteria  for  suitable  aerodynamic  object  functions  can  be  found  In  Ref.  29. 
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Figure  30,  tha  second  airfoil  in  particular,  alao  llluatracaa  another  potential  problem  of  direct 
(inviscld)  wave  drag  minimisation.  In  tha  absanca  of  (direct)  control  over  the  praaaure  distribution  the 
solution  nay  acquire  unrealistically  high  pressure  gradients,  such  as  near  tha  upper  surface  trailing  edge. 

A  strong  point  of  tbs  numerical  optimisation  approach  is  the  possibility  of  selecting  object  functions 
and  constraints  suitable  for  multipoint  designs.  An  example  of  a  two-point  design  problem  directed  towards 
the  design  of  airfoils  with  low  drag  creep  can  be  found  in  Ref.  23.  Low  speed  airfoil  design  applications 
are  considered  in  Refs.  27,  30.  It  is  also  entirely  possible  to  consider,  e.g.  transonic  drag  minimisation 
and  low-speed  stall  requirements  simultaneously. 

While  the  direct  minimisation  of  drag  is  feasible  in  two  dimensions,  it  is  hardly  so,  at  present,  in 
the  case  of  three-dimensional  wings.  Several  unsuecesful  attempts  in  this  direction  can  be  found  in  the 
literature,  (Refs.  24,  25,  26).  The  main  reason  for  this  failure  la  the  lack  of  accuracy  in  the  determina¬ 
tion  of  the  drag  with  the  currently  available  3D  codes  and  the  limited  number  of  mesh  points  (Chapter  4). 
Another  problem  would  seem  to  be  that  the  problem  of  uniqueness  in  three  dimensions  is  even  more  severe 
than  in  two  dimensions.  The  accuracy  problem  nay  be  overcome  when  more  efficient  algorithms  and/or  more 
computer  power  (vector/parallel  machines)  allows  the  number  of  mesh  points  to  be  increased.  The  uniqueness 
problem  would  probably  require  the  introduction  of  more  constraints  or  more  sophisticated  object  functions. 

Summarizing  the  discussion  on  direct  numerical  optimisation,  it  may  be  said  that  the  potential  possi¬ 
bilities  of  the  approach  are  enormous  with,  at  present,  unique  capabilities  such  as  multi-point  and  con¬ 
strained  design.  However,  the  approach  is  alao  unique  in  terms  of  required  computer  resources.  Substantial 
Improvements  in  both  flow  optimization  code  algorlthsui  and/or  computer  efficiency,  relative  to  current 
general  standards,  are  required  before  numerical  optimisation  in  3D  wing  design  can  be  used  on  a  routine 
bases. 

An  alternative  possibility  for  computational  drag  minimization  is  to  use  aerodynamic  (load  and 
pressure  distribution)  rather  chan  geometric  shape  functions  as  independent  variables.  In  this  approach  an 
optimization  algorithm  is  used  to  optimize  Che  pressure  distribution,  e.g.  with  the  objective  to  minimize 
the  drag.  Using  the  latest  available  estimate  of  the  geometry  this  can  be  done  relatively  cheap  through  an 
Induced  drag  (Treffts  plane)  and  a  boundary  layer  code.  While  this  approach  implies  direct  control  over 
Induced  drag  and  viscous  drag  it  does  not  so  with  respect  to  wave  drag.  However,  wave  drag  nay  be 
controlled  indirectly  by  imposing  constraints  on  the  minimum  pressure  or  maximum  upper  surface  Mach  number. 
With  the  target  C  -distribution  established  the  new  geometry  can  be  determined  by  means  of  an  inverse  code. 
Subsequently  the  Bff-deslgn  characteristics  can  be  determined  by  means  of  an  analysis  code.  The  process  is 
repeated  when  the  new  geometry  differs  significantly  from  the  previous  one  or  when  a  geometry  or  off-design 
constraint  is  met.  In  the  latter  cases  (new)  constraints  will  have  to  be  imposed  on  the  values  of  the 
parameters  describing  the  pressure  distribution.  A  flow  chart  of  the  procedure,  which  is  called  inverse 
numerical  optimization.  (Ref.  31),  is  given  by  figure  31. 

An  example  of  a  method  which  fits  into  the  framework  of  Inverse  numerical  optimization  is  one  for 
determining  the  spanvlee  loadings  for  minimum  Induced  drag  on  lifting  parte  of  multiple,  non-planar  config¬ 
urations  (Ref.  32).  The  program  utilises  constrained  optimisation  techniques  and  panel  method  "technology** 
in  the  Trefftz-plane.  It  is  capable  of  dealing  with  pitching  moment  (trimmed  situation)  and  bending  moment 
constraints. 

The  input  required  by  the  program  is  summarized  in  figure  32.  It  comprises  general  geometric  config¬ 
uration  data,  position  of  the  center  of  gravity,  required  (trimmed)  lift  coefficient,  wing,  body  and  tall 
pitching  moment  and  drag  data. 

A  typical  example  of  the  capabilities  of  the  program  is  presented  in  figure  33.  Figure  33a  Illustrates 
the  type  of  configuration  considered,  with  various  vertical  positions  of  the  horizontal  tall. 

Optimum  spanwise  load  distributions  for  a  high  and  a  low  positioned  horizontal  tall  for  a  given  wing- 
body  pitching  moment  and  30  X  MAC  e.g.  location  are  compared  in  figure  33b.  It  illustrates  the  point  that 
optimal  wing  design  cannot  be  achieved  without  taking  account  of  the  tail  configuration. 

Minimum  induced  drag  span  loads  for  two  different  levels  of  wing-body  pitching  moment  for  a  low  tall 
configuration  are  compared  in  figure  33c.  Figures  35d  and  e  illustrate  the  point  that  the  minimum  trlwd 
Induced  drag  of  high  tall  configurations  is  much  more  sensitive  to  variations  in  wing-body  pitching  moment 
and  e.g.  positions  than  that  of  low  tall  configurations.  This,  however,  does  not  necessarily  imply  that  a 
low  tall  configuration  is  to  be  preferred;  there  nay  be  other  design  criteria  that  dominate  this  choice. 

In  conclusion,  figure  33f  presents  the  optimum  span  load  for  a  configuration  which  utilizes  "active** 
loading  of  nacelle/pylon  and  flap  rail  fairings.  In  the  particular  example  a  1.5  X  reduction  of  the  minimum 
trimmed  induced  drag  can  be  realised  by  "active"  rather  than  "passive"  (streamline)  shaping  of  pylons  and 
flap  rail  fairings. 

The  program  has  been  extended  recently  with  an  airfoil  characteristics  data  base  allowing  e.g.  the 
optimization  of  the  spanwise  distribution  of  wing  section  lift  and  pitching  moment  for  minimum  Induced  plus 
airfoil  viscous  drag  for  a  given  spanwise  distribution  of  airfoil  thickness. 

Programs  such  as  just  described  can  be  very  helpful  in  the  conceptual  and  preliminary  design  of  mini¬ 
mum  drag  aircraft.  They  are  also  used  for  creating  the  starting  point  for  the  specification  of  "target" 
pressure  distributions  in  the  detailed  aerodynamic  design  of  wings  by  means  of  Inverse  methods. 

Chordwlse  pressure  distributions  leading*  approximately,  to  minimum  viscous  drag  can.  In  principle,  be 
determined  by  combining  an  optimisation  algorithm  with  a  boundary  layer  code.  Possibilities  along  these 
lines  are  currently  being  explored  at  NLR  using  a  program  system  called  CADOS  (Ref.  33).  The  program  uses 
beslc  "aerodynamic"  as  well  as  spline-type  pressure  distribution  shape  functions  (Fig.  34)  and  can  handle 
laminar  as  well  as  turbulent  flows.  It  can  address  problems  like  the  determination  of  the  pressure  distri¬ 
bution  that,  subject  to  constraints  on,  e.g.  lift,  pitching  moment  coefficient  and  average  pressure  level 
(airfoil  thickness)  leads  to  minimum  viscous  dreg. 

A  simple,  exploratory  example  of  application  is  presented  in  figure  35.  The  figure  shows  the  results 
of  pressure  distribution  modifications  to  a  basic  prsssura  distribution  oy  means  of  spline-type  shape  func¬ 
tions  with  the  objective  to  minimize  the  dreg.  The  boundary  layer  was  not  allowed  to  separate.  Although  the 
problem,  as  formulated,  is  not  necessarily  a  sensible  one,  it  Illustrates  the  basic  idea  of  the  approach. 

In  the  particular  case  the  answer  is,  obviously,  determined  to  a  largo  extent  by  the  constraints  iiqtlled  by 
the  basic  and  spllna  shape  functions.  Other  choices  for  these  will  lead  to  a  different  answer.  The  approach 
is  currently  being  explored  further. 
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It  1*  the  author's  opinion  that  tha  Invars*  numerical  optimization  approach  la  worthy  of  further 
Investigation;  not  In  tha  least  because  the  computational  effort  required  la  much  smaller  than  that  of 
direct  numerical  optimization. 


6  CONCLUDING  REMARKS 


An  overview  has  been  presented  of  the  possibilities  and  problems  associated  with  the  prediction, 
anslysls  and  minimisation  of  drag  using  CFD  methods.  The  discussion  Is  limited  to  sub/transonlc  Invlscld 
(potential  and  Euler)  flow  models,  boundary  layer  and  wake  methods  and  methods  employing  vlacous-lnvlscld 
Interactions. 

From  the  discussion  it  is  concluded  that  the  absolute  ycuracy  of  tha  coyutational  determination  of 
subsonic  aircraft  (wing)  drag  is,  currently.  oT  the  order  of  10  to  30  "counts"  at  best.  The  main  problems 
seem  to  be 

-  lack  of  resolution  due  to  the  limited  number  of  mash-polnts  and  lack  of  convergence  due  to  the  limited 
number  of  iterations  or  time  steps  that  can  be  run  on  the  computer  systems  that  are  currently  available 
In  most  engineering  environments  (particular  Important  for  wave  drag  and.  but  possibly  to  a  lesser 
extent,  for  Induced  drag) 

-  spurious  entropy  production  In  Euler  codes  as  a  result  of  Insufficiently  accurate  Implementation  of 
boundary  conditions 

-  Isck  of  (sppropriate)  methods  and  turbulence  models  for  the  computation  of  wakes  in  pressure  gradient  and 
boundary  layers  with  complex  Interaction  phenomena. 

An  additional  problem  of  the  majority  of  the  available  Invlscld  flow  codes  Is  that  they  provide  values 
for  Integrated  pressure  drag  only  so  that  a  division  Into  wave  drag  and  Induced  drag  Is  not  readily  avail¬ 
able.  In  spite  of  this  situation  It  has  proved  possible  to  devise  engineering  procedures  for  determining 
drag  Increments  for  (closely)  similar  flow  conditions  or  configurations  with  usable  accuracy. 

Closely  related  to  the  problems  mentioned  above  Is  the  fact  that  configuration  optimization  for  drag 
reduction  by  means  of  direct  numerical  optimization  techniques  Is  hardly  feasible  In  current  engineering 
environments,  at  least  for  three-dimensional  problems.  An  alternative  might  be  to  use  numerical  optimiza¬ 
tion  for  designing  minimum  drag  "target"  pressure  distributions  and  to  use  Inverse  methods  to  find  the 
corresponding  geometry  (Inverse  numerical  optimization). 

Summarizing,  It  may  be  stated  the  accurate  numerical  determination  and  minimisation  of  drag  is  current¬ 
ly  almost  "mission  impossible",  at  least  for  three-dimensional  flow.  However,  there  are  exciting  prospects 
for  the  future,  not  In  the  least  because  of  the  continuing  rapid  advancements  In  computer  technology. 
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•  NORMAL  (  PRESSURE  )  FORCES 
— ►  PRESSURE  DRAG 


•  TANGENTIAL  (  FRICTION  )  FORCES 
—•>  FRICTION  DRAG 


Pig.  1  Distinction  of  drag  componanta  according 
to  tppa  of  fore# 


a  FRICTION 

— ►  “VISCOUS”  OR 
(  FRICTION  + 

a  VORTEX  SHEDDING 

— -►  “INDUCED*  OR 


“BOUNDARY  LAYER'  DRAG 
PRESSURE  FORCES  ) 


“VORTEX”  DRAG 


(  PRESSURE  DRAG  ONLY  ) 


I  SUBSONIC  I  DRA6-DUE-T0-L1FT 


CONTAINS: 

.  (  HOST  OF  )  INDUCED  DRAG 
a  (  SOHE  )  VISCOUS  DRAG 
(  a  POSSIBLY  SOME  NAVE  DRAG  ) 


a  USUALLY  ESTIMATED  OR  FITTED  BY  MEANS 
OF  CLASSICAL  INDUCED  DRAG  FORMULA 


Fig.  4  Definition  of  "drag  dun  to  lif  " 


•  SHOCK  WAVE  FORMATION 

— ►  WAVE  DRAG 

(  PRESSURE  DRAG  ONLY  ) 

Fig.  2  Distinction  of  drag  coaponsnts  according 
to  the  responsible  physical  nechaniana 

(  SUBSONIC  )  MINIMUM  DRAG 

CONTAINS: 

•  (  HOST  OF  )  VISCOUS  DRA6 

•  (  SOME  )  INDUCED  DRAG 


•  USUALLY  ESTIMATED  BY  MEANS  OF  FORM  FACTOR  METHODS 


cn  -  c.d  at/ct  ...)  .jss. 

min  F  *r»f 

Pig.  3  Definition  of  "(aubaonic)  alnlaua  drag" 


COMPRESSIBILITY  DRAG 

CONTAINS: 

•  HOST  IF  NOT  ALL  WAVE  DRAG 

•  SOME  VISCOUS  DRAG 

•  SOME  INDUCED  DRA6 

(  DUE  TO  CHANGE  IN  SPAN  L0ADIN6  WITH  MACH  NUMBER  ) 


•  NO  GENERAL  THEORETICAL  BASIS  FOR  PREDICTION 
Pig.  5  Daflnltlon  of  "coapreaalblllty  drag" 
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Fig.  6  Hierarchy  of  nathenatlcal/phyalcal  CFD 
flow  aodeli 


•  (  REYNOLDS  AVE RASED  )  NAVIER  -  STOKES  METHODS 

•  CAN  PROVIDE  FRICTION  DRA6  AND  PRESSURE  DRA6 

•  CONTAIN  AIL  INFORMATION  NECESSARY  TO  DISTINGUISH 
BETWEEN  VISCOUS  DRAG/  INDUCED  DRA6  AND  HAVE  DRAG 

•  GENERALLY  NOT  YET  AVAILABLE  IN  INDUSTRIAL  ENVIRONMENT 

•  INVISCID  METHODS 

•  MODEL  PRESSURE  DRA6  ONLY 


INDUCED 

DRAG 

HAVE 

DRAG 

EULER 

X 

X 

FULL  POTENTIAL 

X 

x*> 

LINEARIZED  POTENTIAL 

X 

- 

•)  APPROXIMATIVE/  ACCEPTABLY  ACCURATE  FOR  Mn<  1.3  ? 

•  BOUNDARY  LAYER  AND  MAKE  METHODS 

•  FRICTION  DRA6  (  BOUNDARY  LAYER  ) 

•  TOTAL  VISCOUS  DRA6  <  MAKE  ) 

Fig.  7  Drag  (component)  prediction  capability 
potential  of  various  CFD  models 


"Sh~'-3 


Fig.  S  Comparison  of  Suler  (FIX) 5 28)  and  full  potential  (FLO 36)  solutions  (Ref.  1) 
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Pig.  9  Example  of  pressure  drag  dependence  on 
■esh  density  (Euler  code,  RAE  2822  air¬ 
foil,  H  *  0.75,  a  °  3°,  strong  shock) 


Pig.  10  Bxsnple  of  pressure  drag  dependence  on 
■eah  density  and  aesb  distribution 
(shock-free  supercritical  airfoil , 

NLR  7301) 


aesb  extent  (Euler  code,  RAE  2822, 
M  *  0.7S,  a  ■  3°,  strong  shock) 


rig.  12  Example  of  pressure  drag  as  a  function 
of  aesb  density  for  a  transport  type 
elag  with  transonic  floe  ((I)FU>22 
prograa) 


Pig.  13  Bxanple  of  akin  friction  prediction 
(I)  (from  ref.  4) 
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Fig.  14  Example  of  friction  drag  prediction 
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Fig.  23  Comparison  of  experimental  (wake  rake) 

and  calculated  viscoua  drag  distribution 
for  transport  type  wing 


Fig.  24  Effect  of  turbulence  model  on  llft/drag 
in  case  of  wake  with  strong  longitudinal 
pressure  gradient 


Fig.  26  Chart  of  relations  between  drag  components 
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Fig.  26  Schematic  of  NLR  XFL022  system  for  the  analysis  of  wing-fuselage  configurations  in  transonic  flow 


b  )  DRAG  CORRELATION  AT  CONSTANT  LIFT 


Fig.  27  NLR  XFL022  drag  analyses  to  narrow-body  transport  configuration 
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Pig.  28  NLR  XPL022  wing- fuselage  transonic  analysis 
system:  examples  of  standard  graphics  output 


Fig.  29  Flow  chart  of  direct  numerical  optlml- 
sat Ion 
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Example  3A-drag  minimization,  M-0.75,  Example  38-drag  minimization,  M-0.75, 

0°  ,  beginning  analysii  No.  27.  cu 0°  ,  beginning  with  analyait  No.  48. 


Fig.  30  Example  of  non-uniqueness  of  wave  drag  minimisation  problem 
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